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Heat dissipation in current-carrying cryogenic nanostructures is problematic because the phonon 
density of states decreases strongly as energy decreases. We show that the Coulomb interaction 
can prove a valuable resource for carrier cooling via coupling to a nearby, cold electron reservoir. 
Specifically, we consider the geometry of an electron bilayer in a silicon-based heterostructure, and 
analyze the power transfer. We show that across a range of temperatures, separations, and sheet 
densities, the electron-electron interaction dominates the phonon heat-dissipation modes as the main 
cooling mechanism. Coulomb cooling is most effective at low densities, when phonon cooling is least 
effective in silicon, making it especially relevant for experiments attempting to perform coherent 
manipulations of single spins. 



I. INTRODUCTION 

As researchers continue to probe smaller electronic de- 
vices at lower temperatures, a detailed understanding 
of heat management applicable on such length and en- 
ergy scales becomes increasingly important. For exam- 
ple, recent experiments to detect the spin resonance of 
a single electrorP and to perform fast charge sensing in 
few-electron quantum dotfl are both limited by heating 
effects. Other applications, such as the search for the 
v = 5/2 non-abelian quantum Hall state, are expected to 
require very low temperatures,'^ making the development 
of schemes for cooling such devices a necessary challenge. 

The main problem is that whenever current is applied 
to a device to perform a measurement, the conduction 
electron temperature increases due to Joule heating.l^ In 
devices operating near room temperature, heat can be 
readily dissipated through phonons, as the conduction 
electrons and lattice are strongly coupled. However, as 
temperature is decreased, the conduction electrons de- 
couple from the lattice. The phonon modes contribute 
less and less to cooling because the phonon density of 
states decreases as energy is decreased.^l Hence, as the 
system gets colder, it becomes more difficult to cool via 
conventional means. 

Besides phonon cooling, systems can be cooled by elec- 
tron diffusion through the leads.'^ However, in common 
nanoscale devices, the leads extend hundreds of microns 
from critical regions to regions that are well cooled. This 
large distance scale limits the effectiveness of electron 
diffusion for device cooling.'Sl 

Here, we investigate using the Coulomb interaction di- 
rectly, a strategy for cooling that remains largely un- 
explored. That is, we consider placing a cold conduc- 
tor nearby the hot conduction electrons, and using the 
Coulomb interaction for heat transfer. The cold body 
could then be heat-sunk directly, without having to worry 
about interfering with the operation of a device. 

While it may seem that remote Coulomb interactions 
are not strong enough to facilitate meaningful power 
transfer, several recent experiments have shown that re- 
mote interactions can indeed drastically affect electron 



relaxation. For instance, the widely studied Coulomb 
drag (CD) effectPHSl involves the transfer of momentum 
from one two-dimensional electron gas (2DEG) to an- 
other via the Coulomb interaction, due to the layers' 
close proximity. Another example of the importance of 
remote Coulomb interactions arises in the metal-oxide- 
semiconductor (MOS) geometry, where it has been found 
that device performance can be reduced due to interac- 
tions of the conduction electrons with those in the gate 
when the distances are too small.ti^E^ 

In this paper, we consider two parallel, silicon 2DEGs, 
separated by tens of nanometers. We make the simplify- 
ing approximation that both 2DEGs are of zero thickness. 
Similar devices have been implemented e xperi mentally in 
the form of silicon electron-hole bilayers .^^HAI One of the 
layers is taken to have a temperature on the order of tens 
to hundreds of mK; this is the active layer that we are 
interested in cooling. The second layer is taken to be a 
heat sink, and is assumed to be at temperature T = 
K. Although cooling this heat sink layer would be a chal- 
lenging engineering problem, electron diffusion through 
grounded, close proximity leads can effectively cool sam- 
ples to about 10 mK.'-^^ Since the two layers are electri- 
cally decoupled, this would not interfere with electrical 
measurements on the active layer. 

We study the temperature, separation, and carrier den- 
sity dependencies of the heat transfer, and compare it to 
experimental results for heat dissipation due to phonons. 
The Coulomb interaction is found to be competitive and 
even dominant over phonons for a range of tempera- 
tures and densities, for separations up to several tens 
of nanometers. Specifically, we find that lowering the 
electron density enhances the power transfer due to the 
Coulomb interaction, but decreases the power dissipation 
due to phonons. 

This paper is organized as follows. Section [H] uses a 
formalism similar to that used for CD to formulate the 
physical problem. Next, Sec. |III| finds an expression for 
the power transfer and discusses its asymptotic behavior 
in both near and far distance regimes. 
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pares the results of the previous sections to experimen- 
tal results for phonon-mediated cooling. Finally, Sec. |V] 
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discusses the implications of cooling nanosystems using 
carrier-carrier interactions, and suggests directions for fu- 
ture study. 



II. FORMULATION OF THE PROBLEM 

In this section, we describe the physical situation we 
will consider throughout this paper: two parallel 2DEGs, 
one serving as a heat sink for the other. We then review 
the standard scattering formalism that is used to per- 
form Coulomb drag (CD) calculations. Using this Boltz- 
mann transport formalism, we write down an equation 
for power transfer, which we will evaluate in subsequent 
sections. 

The physical situation we consider here is very sim- 
ilar to th at w hich has been well-studied in the CD 
literature .EtSllSI Specifically, we have a sample that con- 
tains two 2DEGs that are parallel but spatially sepa- 
rated by some distance d. In the case of CD, one of the 
2DEGs is driven by a current, while the other is not. Be- 
cause of the current flow, the distribution function in the 
first layer is out of equilibrium, and the resulting charge 
fluctuations generate a response in the second layer, me- 
diated by the Coulomb interaction between the layers. 
Although complicated by screening, the basic picture is 
that the presence of the current in the first layer "drags" 
electrons in the second layer, creating a net voltage. 

In our case, we do not consider the non-equilibrium ef- 
fects of a current flowing. Rather, we suppose that each 
layer is internally at thermal equilibrium, but at differ- 
ent temperatures. We will consider one layer (the active 
layer) to be at a finite temperature T, and the second 
layer (the heat sink) to be at T = K. Intuitively, we 
expect that energy should be transferred from the active 
layer to the heat sink. Microscopically, this is due to den- 
sity fluctuations in the hot layer causing responses in the 
cold layer, mediated by the Coulomb interaction. For- 
tunately, we can borrow much of the initial setup of the 
problem from the CD formalism. However, the evalua- 
tion of the resulting expression is quite different because 
we focus on energy between 2DEGS of unequal temper- 
ature, rather than momentum transfer due to a driving 
electric field. 

Since we will be working at very low temperatures (< 1 
K), binary Coulomb collisions are the primary method 
of heat transfer between layers. If we were to consider 
temperatures T > 0.2Ep/kB, where Ep is the Fermi en- 
ergy, we would also need to take into account collective 
scattering effects, the so-called plasmon enhancement .^^^ 
Hence, we consider interactions that transfer energy from 
an electron in the active 2DEG (layer 1) to an electron in 
the heat sink 2DEG (layer 2). The electrons involved in 
the interaction have initial (2D) momenta (ki, k2) and fi- 
nal momenta k2) = (ki -I- k, k2 — k), where k is the 
transferred momentum. This carrier-carrier scattering 
falls into the category of distinguishable particle scatter- 
ing since events only occur between particles in different 
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FIG. 1. (Color online) The dimensionless function Y{ri,(o) 
(Eq. ([3|) plotted over rj for various values of and with layer 
densities ni = n2. Here, rj — k/kp is a scaled momentum 
and ^0 = ksT/Ep is the characteristic magnitude of energy 
fluctuations due to temperature. The function follows a power 
law for both low and high r/, with a maximum occurring for 
an intermediate rj, which we denote as ry*. 



layers. The formalism for treating such a scattering prob- 
lem is well known,!^ and the power transfer is shown in 
appendix [K\ to be 



16^3 
(2^ 



d^ki(fk2(fk ■ E-T, 



(1) 



where A is the sample area, E is the transfer energy for an 
individual event, and F is the scattering rate. In the next 
section, we evaluate Eq. ([T]) using the Thomas-Fermi ap- 
proximation to describe screening (which in this context 
is equivalent to the RPA approach taken in Ref . [TB]) . 



III. RESULTS 

In this section, we evaluate Eq. ([T]), using some stan- 
dard methods incorporated into the calculation of CD. 
However, the resulting expression is quite different, so 
we work through its derivation in appendix [C] We set 
the temperature of the active layer to T, the tempera- 
ture of the heat sink layer to absolute zero, the Fermi 
level of the active layer to Ep, and the Fermi level of the 
heat sink layer to Ep/x, where x = ni/n2, the ratio of 
carrier densities between layer one and layer two. The 
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power transfer is then 
where 



J] 



smhiri/Tjo] 



YivXo), (2) 



Y{vXo)= I dC^[coth(C/Co)-l] 



X Re 



(v/2(2 + C)772-774_^2 




(3) 



and the dimensionless parameters are C = E/Ep, 
T] = k/kp, Co = ksT/Ep and rjQ = \/{kpd)^ with 
kp = yj2m* Ep /h the Fermi momentum. Here, the 
•q/ sinh(77/?7o) term is due to the intcrlayer Coulomb inter- 
action (Eq. ( |A5| ), while the distribution functions give 
rise to Y{'q, Qq). 

The function Y{rj, Co) is plotted in Fig. [l] where it is 
shown that K is a peaked function in rj varying as a power 
of rj on either side of the peak. We define 77* to be the 
location of the peak, and note that 77* ~ Co- Physically, 
y(?7, Co) tracks the availability of energy fluctuations cor- 
responding to a particular momentum transfer k = rjkp 
and temperature fc^T = CoEp. If 77 < 77*, K is limited by 
the Fermi-Dirac distributions that govern the occupation 
of states in each 2DEG. For 77 > 77*, F is instead con- 
strained by the temperature difference between 2DEGs. 

From Eq. (|2| , we see that the Coulomb potential causes 
77 to be cut off at approximately 779. Hence, there are two 
asymptotic regions of interest: when 770 <^ rf and when 
Vo ^ V* ■ III the first region, the Coulomb potential cuts 
off the integration over 77 well before rj* , which corre- 
sponds to large separations between 2DEGs. Here, the 
separation distance limits the magnitude of the momen- 
tum transfer, which in turn limits power transfer. For 
% 3> 77*, the Coulomb interaction truncates the 77 inte- 
gration after 77*, corresponding to small separation. In 
this regime, Y is already rapidly decreasing, so the power 
transfer is instead mainly constrained by the temperature 
difference between 2DEGs. The crossover between these 
two regions occurs when 770 « 7/* ~ C07 corresponding to 
a separation of d « Ep / {kpksT). 

In appendix [D] we work out the asymptotic forms for 
power transfer. We consider the specific case of equal 
density 2DEGs, when x — 1. In the large separation 
regime when d S> Ep / {kpksT), we find 



P 



kpti^ 



keT 
d^ 



[8.3 + 13.0- log (fc^d)], (4) 



where we use (~) to denote asymptotic equivalence. For 
the short distance limit where d <C Ep/{kpkBT), we 
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FIG. 2. (Color online) Calculated values for the power per 
unit area P/A transferred via the Coulomb interaction ver- 
sus separation (Eq. ([2])). Here, the sheet density is 4 x 10^^ 
cm~^ in both layers, the temperature of the active layer is 
T = 100 mK, and the temperature of the heat sink layer is 
K. As can be seen by the dashed lines, the power trans- 
fer varies as approximately at small distances and 
at large distances. The crossover length scale occurs when 
d ^ Ep/ikFkBT). 
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l2?,Epm* V 
0.46 - 1.32 • log 
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knT 



Ep 



0.81 - log (kpd) 



(5) 



Hence, up to logarithmic corrections, P/A oc T/d^ for 
large distances and P/A oc /d^ for small distances, 
which can be qualitatively understood as follows. At 
low temperatures (T ^ ksEp), it is reasonable to as- 
sume that energy fluctuations are small and concentrated 
about the Fermi level, so that the transfer momentum 
obeys k <^ kp. Expanding Eq. ([s]) for small 7/, while 
working in the large separation regime where rj <^ rj* , we 
find that 



(6) 



Here, the scaling is determined by the Fermi-Dirac distri- 
butions limiting the power transfer. Likewise, if we work 
in the small separation region where 77 3> 77* , wc find 



Yiv,Co)'x 



CI 



(7) 



4 



where the scahng is now determined by the layer tem- 
perature. In these hmits, since the Coulomb interaction 
sets the scale of 77 oc 770, we can easily see the rough 
dependences (neglecting the logarithmic corrections) via 
power counting in Eq. ^ . Fig. [2] shows the numerical 
evaluation of Eq. ([2]) as a function of separation, clearly 
demonstrating both distance regimes. 



IV. COMPARISON WITH COOLING DUE TO 
PHONONS. 

In this section, we compare Coulomb-mediated cool- 
ing to experimentally measured energy dissipation due 
to phonons at low temperatures. In silicon-based het- 
erostructures at low temperatures, two types of phonon 
couplings are important: acoustic phonons governed 
by a deformation potential coupling,'^^ and the Pekar 
coupling .fi^ Pekar phonons arise from the sharp electro- 
static confinement potentials present in heterostructure 
devices, such as quantum wells, and hence are not present 
in bulk samples.^- They also share a characteristic de- 
pendence with piezoelectric phonons,'^ making them es- 
pecially important in low-temperature experiments with 
non-polar materials, such as few-electron quantum dots 
in Si. 

Indeed, both deformation potential and Pekar phonons 
have been experimentally observed in silicon-based het- 
erostructures at low temperature.'^ The characteristic 
temperature dependence for deformation potential cou- 
pling is T^^^ so for very low temperatures we expect 
Pekar phonons to dominate, while for higher tempera- 
tures deformation potential phonons become more im- 
portant. 

As established in Eq. ([5]), for small separations the 
power transfer to the heat sink layer via the Coulomb 
interaction varies as T^. Whether or not this Coulomb 
cooling is larger than phonon cooling over a given tem- 
perature range depends on the numerical magnitude of 
Eq. ([2]), which we now compute. We compare Coulomb 
cooling to experimental measurements of phonon me- 
diated cooling in Ref. [HI There, it is found that the 
power dissipation due to phonons is Pph/A — aT^ + bT^, 
where a = 2.2 x 10"® W K'^ cm'^ and 5 = 5.1 x 10"^ 
W cm^^. The structure used is a silicon MOS inver- 
sion layer, with dielectric thickness 200 nm and carrier 
density 5.4 x 10^^ cm~^. 

It is known that the phonon couplings depend on the 
electron density, with P cx n^'^/^ for the deformation 
potential coupling.^ Pekar phonons have both an explicit 
n~^/^ dependence and a dependence on the electric field 
at the 2DEG of ^^^.EIl since in a 2DEG F cx Pekar 
phonons scale as P cx n'^/^ in total. By comparison, 
Eq. ([5| tells us that for equal density 2DEGs, Coulomb 
power transfer goes like P oc 1/n. 

Typically, in an experiment the density is fixed by 
desired electronic properties (for instance, the ability 
to pinch off current with depletion gates). For low- 
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FIG. 3. (Color online) Calculated values for the power per 
unit area P/A transferred between an active (T = 50 mK) 
and heat sink (T = K) 2DEG via the Coulomb interaction 
as a function of the layer density at three different values for 
the separation between layers. Here, the densities of both 
layers are identical. For comparison, P/A due to phonons 
from experimental data in Ref. [G] scaled for changing density, 
is shown as a dotted line. 



temperature applications that attempt to reach few- 
electron regimes, it is desirable to have a low density. 
It is therefore important to determine the dependence of 
power transfer on the density n. Fig. [3] shows the effect 
of varying the layer density on the power transfer, for 
three different layer separations at constant temperature 
T = 50 mK. As expected, the Coulomb power transfer is 
greatest in the case of small density, making it especially 
pertinent for few-electron experiments. 

It is important to note that our formalism for static 
screening is only valid when the transfer momentum 
obeys k < 2^^- ,1221 which means that our we cannot make 
the density too small. The Coulomb interaction limits 
the transfer momentum to fc < l/d. Hence, setting k = 
1/d for d = 10 nm corresponds to n > 0.79 x 10^^ cm^^. 
Another constraint on low-density 2DEGs is the metal- 
insulator transition, which occurs for sufficiently low den- 
sities. In silicon MOS structures, the critical value of 
density is known to be around Uc ~ 1 x 10^^ cm~^.l^ 
More recently, calculations for dopantless Si/SiGe de- 
vices predict that this value can be much lower, about 
« 2 X IQio cm-2.22 

In Fig. |4] we plot the temperature dependence of the 
power transfer per unit area P/A for several separations, 
and compare with the power dissipation due to phonons. 
There, we fix the carrier density to be n = 1 x 10^^ cm"^ 
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FIG. 4. (Color online) Scaled values for the power per unit 
area P/(Ar*) transferred via the Coulomb interaction be- 
tween two 2DEGs versus temperature at three different val- 
ues for the separation between the layers. For comparison, 
the scaled P/ (AT*) due to phonons from experimental data 
in Ref.|B]is shown. The density of both layers is 1 x W^^ /cm^ . 



for both layers. One sees that for small separations (less 
than 20 nm), Coulomb- mediated power transfer exceeds 
phonon power dissipation over a potentially wide tem- 
perature range (roughly up to 300 mK for a 10 nm sep- 
aration) . 



netic formulation by Volokitin and Persson.l^^I^Il How- 
ever, these are largely interested in hot devices, where 
complicating features such as plasma excitations are im- 
portant. Further, as noted in Ref. there are dis- 
crepancies between this electromagnetic formalism and 
Boltzmann transport approaches. More recent work by 
Kriiger, Emig, and Kardar extends the electromagnetic 
formalism to arbitrary geometries with a focus on heat 
transfer.l^fil It would be beneficial to compare the present 
work to the electromagnetic treatments to attempt to 
address the origin of any discrepancies. 

While the results for two parallel 2DEGs are promising, 
one could almost certainly engineer a better geometry for 
optimizing heat dissipation. Indeed, the main reason for 
a preliminary evaluation of the 2DEG-2DEG geometry 
was due to its computational simplicity. An idea for a 
more effective heat sink might be a standard MOS geom- 
etry, or a top-gated nanostructure. Due to the drastically 
higher density of states in the metal, one could expect an 
enhanced power transfer. However, screening would also 
be enhanced, so careful calculations, similar to those pre- 
sented in this paper, should be done for that geometry. 
Also, studying the effects of high-k dielectrics might be 
fruitful, since the power transfer scales as e^. 
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Appendix A: Derivation of the power transfer rate 



V. DISCUSSION 

Understanding relevant heat dissipation mechanisms 
at low temperatures in electronic devices is an important 
problem, especially as spin-based, few-electron devices 
mature. In this paper, we considered a geometry con- 
sisting of parallel 2DEGs in silicon and calculated the 
expression for power transfer between two layers at tem- 
peratures T > and T = respectively, in the approx- 
imation of Thomas-Fermi screening. We then presented 
analytic results for the asymptotic regimes of small and 
large separations. We showed that in this geometry, 
power transfer due to the remote Coulomb interaction 
can be the dominant heat loss mechanism. This Coulomb 
cooling is most effective at low densities, making it es- 
pecially important for experiments attempting to access 
few-electron regimes. 

There have been a number of studies of heat transfer 
between close bodies, including a semiclassical kinetic 
treatment by Boiko and Sirenko^^ and an electromag- 



In this appendix, we briefly sketch the derivation of Eq. 
([T]), the formal expression for the power transfer between 
two 2DEGs, using the methods of Ref. [T71 First, recall 
that we are interested in the scattering of two particles 
with initial (2D) momenta (ki,k2) and final momenta 
(k']^,k2) — (ki+k, k2 — k). The transition rate F for 
the above process is given by the balance equation 

r (ki, ka; k'l, k^) ^ S (ki, ka; k'l, k^) (Al) 



f(2) 



where S is the transition rate given that the appropriate 
states are available, f^-^^ is the Fermi-Dirac distribution 
function in layer one. 



1 + exp 



Ek — Ef 

knT 



(A2) 
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where i?k = h^k"^ / {2m*), Ep is the Fermi level, m* is 
the effective mass, and T is the temperature in layer one, 
and /j^ is likewise the Fermi-Dirac distribution function 
in layer two. Note that here we restrict our attention to 
Fermi-Dirac distribution functions, but Eq. ( Al ) remains 



valid even for non-equilibrium distribution functions. 
The first term in the square brackets of Eq. ( Al I can 



be understood as the particles starting with momenta 
(ki,k2) and ending with momenta (k']^,k2). The second 
term corresponds to scattering from momenta (k'j^, kg) to 
momenta (ki, k2). To calculate the scattering rate 5, we 
use Fermi's golden rule: 



5(ki,k2;k;,k'2; 



27r 



\H\ 



(A3) 



where H is the interaction matrix element, 

= |^.cr(fc,d), (A4) 

and (pscr is the Fourier transformed screened Coulomb 
interaction between layers. Defining the Thomas-Fcrmi 
screening wavevector kxr = ^.m* /{■kH'^) ■ g'^/(2eoeh), the 
Fourier-transformed screened Coulomb interaction be- 
tween the layers, calculated within the Thomas-Fermi 
approximation ii 



4>scrik,d) 



kq 



Ak^p£Qei, smh(kd) ' 



(A5) 



where we have assumed k fcrF- For clarity of pre- 
sentation, we present a self-contained derivation of this 
expression in appendix [B| 

Now that we have an expression for the scattering rate 
between particular states, we obtain the power transfer 
between the layers: 



P^16 J2 £^-r(ki,k2;ki -Hk,k2 -k), 



(A6) 



ki,k2,k 



where the factor of 16 is due to spin degeneracies of two 
and valley degeneracies of two in each electron layer,'^^ 
and E = -Eki+k~£'ki i the transferred energy. Converting 
the sum to an integral gives us Eq. ([I]). 



Appendix B: The screened Coulomb potential 

In this appendix, we present a self-contained deriva- 
tion of the screened interlayer Coulomb potential within 
the Thomas-Fermi approximation. Although this result 
can be obtained as a special case of the random phase 
approximation result as described in Ref. (THl assuming a 
static screening formalism from the beginning results in 
a considerably more transparent calculation. The tech- 
nique we present here can also easily be implemented 
numerically to treat more complex geometries. 



To start, we consider placing an electron into one of the 
2DEGs. This results in an external, unscreened potential 
(j)ext{^,z) due to the external electron, where r is the 
2D position within the plane of the 2DEG. The electron 
gas in both layers can rearrange to screen this external 
charge, resulting in an induced potential cjiindi'*^, z). The 
screened potential that an electron in the other layer feels 
is then 4>scr{Y, d) = ipexti'r, d) + 0md(r, d), where we have 
assumed that the two 2DEGs are separated by a distance 
d. Our objective is to calculate (l)ind, from which we can 
compute (pscr- 

We assume that our system is translationally invariant 
in the plane parallel to the 2DEGs, which we define to 
be the x — y plane. It is convenient to exploit this trans- 
lational invariance by taking a Fourier transform of the 
Poisson equation in the x — y plane, yielding 

{dz€{z)dz - e(z)fc^) 0„id(k, z) = -p„jrf(k, z), (Bl) 

where we denote the Fourier transform of a function 
/(r, z) as 



d\f (r, z) 



-irk 



(B2) 



In Eq. (Bl ), pind is the induced charge density, responsi- 



ble for the production of cpind, and e{z) is the dielectric 
function. In a homogeneous medium, Eq. (Bl) has the 
general solutiorP^ 



(i>ind{k,z) 



1 



2keoeb 



dz'e-'^l--^'lft„rf(fc,z'), (B3) 



so to find (pind, we must calculate pind- 

To determine Pmd, we first note that the total charge 
density ptot obeys 



Ptotir, z) ^ po+ Pind, 



(B4) 



where po is the charge density without an external charge 
present, and we have neglected the small density contri- 
bution from the external charge itself. The dispersion 
relation for the electrons is given approximately by 



2l2 



Eik,r,z) 



ffk 
2m* 



- q(t>scr{-r,z), 



(B5) 



where — g is the charge on an electron. By using the 
functional form of the Fermi-Dirac distribution, we can 
view charge density as a functional of Fermi energyP 

/5md(r, z) « Po {Ep + q(j>scr) - Po (Ef) ■ (B6) 
cr ^ Ep, to first order in 4>scr 



Now, assuming that 
Eq. ([B6| is 



Pind{r,z) « -q- 



driQ 

Ie 



(B7) 



where po = —quQ. For low temperatures, dno/dE^^ « 
g{r,Ep), the local density of states evaluated at the 
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Fermi level, which might vary spatially. For our geom- 
etry with two 2DEGs separated by a distance d, g only 
varies in the z direction: 



g{z) ^ 92D {S{z) + S{z - d)) 



(B8) 



where g2D is the energy-independent two-dimensional 
density of states. 



Substituting Eq. p37|) into Eq. p33|), we find 

(f>indiz) 



2fceoe6 

+ e-^l^-'^l (^4>^nd{d) + 4>eM 



(B9) 



The external potential due to the external electron in the 
first layer satisfie^^Hl 



4'extik, z) 



(BIO) 



Hence, evaluating Eq. (B9| for z = and z = d leaves 



us with a system of two linear equations. Solving gives 
Eq. ( A5 1 , where we note that the Thomas-Fermi screen- 



ing wavevector krF is defined to be twice what is typical 
for GaAs, due to the extra valley degeneracy in Si, and 
we have assumed that k <^ kxp- 



Appendix C: Calculation of the power transfer 

This appendix presents the derivation of Eq. ^ . The 
calculation begins similarly to those done in the case of 
Coulomb drag.fi^ However, it proceeds quite differently 
because the symmetry of the momentum transfer rele- 
vant to CD differs from that of power transfer, which we 
consider here. Following the CD literature,'^^ we seek to 
decouple the fci and fc2 integrals. First, it is conventional 
to split the energy-conserving delta function in Fermi's 
golden rule by introducing an integration over the trans- 
fer energy. The relevant identity is^^ 

S (-Eki + Ei^2 ^ -S'ki+k — -Ek2-k) (CI) 

= J dES{E + Sk, - E^^+],)5{E - Sk, + £^k,-k). 

Next, note that Fermi-Dirac distributions at a common 
temperature T satisfy the algebraic relationship for two 
energies and E^^ 



f{E^){l- f[E^ + Ey)) 



1 - e-B„/(fcBT) 
A third, useful algebraic identity is 
11 11 



1 - e-° 1 - 6+'' 1 - e+° 1 - e-'' 



— coth b — coth ( 



(C3) 



which is verified by using the definition of the hyperbolic 
tangent. With these identities, it is tedious but straight- 
forward to show that Eq. ([T]) can be written as 



P 
~A 



where 



16g2 
coth 



d^kdE ■ E 
E . 

coth 



kBT2 



tot{k,d) I(k,E)J(k,E) 
E 



ksTi 



(C4) 



I{KE)= J d^kid{E + E^,-Ek,+^) 



f^'HE^J-r^>iE^,+E) 



and 



J(k, E)^ J d^aSiE - Ek, + i;k,-k) 



(C5) 



f^'>{E^J-f'>iE,,,-E) 



(2) 



(C6) 



Again, /^^^ is the Fermi-Dirac distribution function of 
layer 1 and /^^' is the Fermi-Dirac distribution func- 
tion of layer 2. We next make the simplifying assump- 
tion that the carriers in the two layers have have the 
same effective masses, but possibly different Fermi lev- 
els. Assuming that the temperature is sufficiently low, 
we also approximate the distribution functions as step 
functions at the Fermi level, from which it follows that 
/(k, E) ss — J(k, E) when the Fermi levels are the iden- 
tical. 

We calculate /(k, E) within the effective mass ap- 
proximation with a simple parabolic dispersion, E — 
h^k"^ / {2m*), by using Cartesian coordinates and integrat- 
ing over ki, yielding: 



HKE) = 



h^k 



(C7) 



where m* is the (transverse) effective mass, Ea — Ep — 
Eq — E, Ep — Ep — Eq, Ep is the Fermi level of layer 1 
and 



E. 



2m* 



Em* 
l\?k 



(C8) 



It is now useful to switch to dimensionless coordinates, 
where we define C, = E /Ep and rj = k/kp, where kp = 
\/2m* Ep/h is the Fermi momentum. Doing this gives 



rRe 



(v/2(2 + c)v^ -v'-e 

- ^2(2 - Orj^ e) ■ (C9) 



The calculation of J is very similar, except that the Fermi 
level of layer two is taken to be Ep/x, where x = n\jn2 
is the ratio of carrier densities. Recall that we wish to 
consider systems where Ti ^ T2. Hence, for simplicity 



8 



we let T2 « 0. Introducing the parameters Co = ksTi/Ep 
and i]o = l/{kFd) and substituting Eqs. (A5) and (C9) 
into Eq. (|C4|) gives us 



A 6Ah \ g2 



drj 



^sinh(77/7?o) 

which is Eq. where Y is defined by Eq. ^ 



Yin, Co 



(CIO) 



Appendix D: Asymptotic analysis of heat transfer 

In this appendix, we seek to obtain an accurate 
analytic expression for Eq. ^ for both large {d 
Ep/ikphBT)) and small {d <C Ep/(kpkBT)) separations 
between the 2DEG layers. In the following, we set the 
densities of the two 2DEGs to be equal for simplicity. 
To, proceed, we first expand for small momentum exci- 
tations about the Fermi level. Using this, we work out the 
asymptotic form of Y (Eq. (|3])) on either side of its peak. 
Then, we calculate the resulting integral in Eq. ([2]), and 
derive formulas for asymptotic power transfer for both 
large (Eq. Q) and small (Eq. ([5])) separations. 



1. Asymptotic forms of Y 

To begin, we calculate the asymptotic forms of Y . 
Since we are at low temperatures, we may assume that 
the transfer momentum k <^ kp, the Fermi momentum. 
In this approximation, we find that 



2"-"' c(coth(c/co) - 1) ^ev 

4 
2"+"' C(coth(C/Co) ^ 1) 



The two limiting cases we consider are when ?7 ^ Coj 
corresponding to the region well to the left of the peak 
in Y, and ?7 S> Coi corresponding to the right of the peak. 
When rj <^ Co, we make the approximation that 



coth(C/Co) « Co - C- 
Using this, Eq. (|D1|) reduces to 



YivXo) - 4Colog(4/r;)r?2. 

'?<Co 

When 77 ^ Coi we may take rj ^ (, so 
In this limit, the second integral in Eq. (Dl) does not 



(D2) 
(D3) 

(D4) 



contribute. Evaluating the first integral, we find 

/■A poo 

YivXo) — / dx(cothx- 
')»Co V Jo 



(D5) 
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FIG. 5. (Color online) The rescaled function v(e) (Eq. (|D7|)), 
plotted versus the scaled coordinate = rj/C,o- Here, 77 = 
k/kp, the momentum transfer scaled by the Fermi momen- 
tum, and Co ~ ksT/ Ef- For small 9, v takes on the expected 
asymptotic value of one. The curves essentially coincide until 
around >^ 1. This enables us to treat v as approximately 
independent of Co before that point. 



which reduces to 



Y[v.Co. 



r,»Co 77 120' 



(D6) 



Now that we have the asymptotic forms of Y on either 
side of the peak, we may proceed to evaluate Eq. ([2| in 
.the limits of 770 ^ Co and ?7o ^ Co- 



2. Power transfer in the limit of large separation 

Next, we evaluate Eq. ([2]) in the large-distance limit, 
when 770 <C Co: which corresponds to c? 3> Ep/{kpkBT). 
To do this, we first define a scaled function v. 



viOXo) 



Y{9Co,Co) 



4C3log(4/(^^Co))^^' 



(D7) 



which is just Y scaled by its asymptotic value in the 
region where 77 ^ Co a function of the scaled coordinate 
9 = r]/(o- We plot v for various values of Co in Fig- [s) 
The integral we need to evaluate can be written as 



A(r;o,Co) = < 



.o'-^^m^MM^ (D8) 



sinh^(0Co/'7o) 



which is related to Eq. [2] by 



J^p I £0^6 \ A / /- \ 

' ' A('7o,Co)- 



P^E^ ( epeb 
A 64/7 V 1^ 



(D9) 



9 



Since we are in the region where 770 Co, we can approx- 
imate v{9,(o) ~ 1. The integration can then be carried 
out numerically, resulting in 



HvoXo) Co^^ (8.3 -13.0 log 7?o), (DIO) 



which reduces to Eq. Q when inserted into Eq. (D9). 



3. Power transfer in the limit of small separation 



Now, we evaluate Eq. ^ in the small-distance limit, 
when 770 » Co, corr esponding to d <^ Ep/{kpkBT). We 
begin with Eq. ( D8 1 , but unlike before we cannot assume 
1. Instead, we note that from Fig. [s] 



that v{9Xo) = ^- instead, we note tnat trom t ig. [Oj v 
is approximately independent of Co until some cutoff 9, 
6'c > 1. Hence, we split the integration region of A into 
two pieces at 9c- For 9 < 9c, we take v to be independent 



of Co: and also 



sinh(6'Co/?7o) 



9_Co 
Vo ' 



(Dll) 



where this second approximation is valid since 9 < 9c <^ 
ryo/Co- For 9 > 9c, we calculate v{9,(q) according to the 
asymptotic formula for Y in the limit of 77 3> Co, given in 
Eq. ( D6 ) . These approximations result in 



A(77o,Co) KoVo / d9 \ogi4/i9Co))9M0) 



120 



So 



d9- 



sinh (6'Co/?7o) 



(D12) 



In the limit where Co ^ ^70 and with 9c = 2.0, the result 
is 

A(%,Co) - Co%' (0.458 - 1.32 log Co + 0.182 log 770). 

(D13) 

In this evaluation, we picked Co = 0.001 for the calcu- 
lation of V in the region where it is approximately Co- 
independent. Inserting this result into Eq. (D9| gives 
Eq. ([5|, as desired. 
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